Computation of pion and kaon heavy ion multiplicities in a gluon-meson model 
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In high energy Heavy Ion Collisions, the onset of the quark-gluon plasma is the colour glass 
condensate, dominated by gluons. The final state is hadronic, and dominated by pions and kaons. 
Here we investigate an effective approach of QCD with these bosonic fields and which can help to 
describe the transition of gluons into light mesons. Formally, our approach consists in integrating 
out the quark fields from the QCD path integral. In this way the fermionic fields are replaced 
by light mesons, such as the pions and sigma field. We apply our effective action to compute the 
number of pions and kaons per gluon emitted by a Boltzmann gluon gas, their multiplicities as a 
function of the gluon mass. We conclude that an effective gluon mass remains finite at T = Tc- 



I. INTRODUCTION 

In this work we develop a Lagrangian where only 
bosons, i. e. gluons and mesons, are the active degrees 
of freedom, and apply it to study the multiplicities in 
heavy ion collisions, the number of pions and kaons per 
gluon, and the gluon mass at the onset of the deconfine- 
mcnt/confinement phase transition 

Effective approaches of QCD have been widely used to 
study the properties of strong interactions 0. Quark 
models, meson effective models, or models combining 
both quarks with mesons are used thoroughly to explore 
hadronic physics. Although gluons have been proposed 
already in the 70's together with the theory of strong 
interactions, QCD, in effective models it is common to 
assume that the gluons are integrated out, and only con- 
tribute indirectly through the quark or effective hadron 
interactions. 

Nevertheless, there are two rapidly developing QCD 
domains where gluons are either easier to work with, or 
are phenomenologically more relevant. In many-body 
systems, the Grassmann variable nature of the quarks 
makes them technically much more difficult to address 
than the bosonic gluons. In particular. Lattice QCD first 
developed and applied pure gauge or quenched techniques 
since working with dynamical quarks is computationally 
very expensive Moreover, in high energy Heavy Ion 

Collisions [J , it was proposed successfully that the onset 
of the quark-gluon plasma is the colour glass condensate, 
dominated by gluons The final state is hadronic, 

and dominated by pions and kaons. For instance, in the 
many particle BAMPS set-up for heavy ion collisions, 
fiot the simulations are performed with gluons only 
form the onset, and mesons arc included as final states 
of the hadronization. From the QCD perspective, in this 
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case, quarks and not gluons are integrated out and effec- 
tive mesons are included. Note also the recent work of 
Weinberg where an effective Lagrangian with gluons, 
in addition to pions and quarks, was put forward. 

Formally, the different effective approaches to QCD 
can be seen as the result of integrating out some degrees 
of freedom from the QCD Lagrangian For instance, 
when only gluons are integrated out, one obtains a NJL- 
like theory [l^l or a quark model [l^. If, in addition, 
also quarks are integrated out from the NJL model, one 
is left with a purely linear a model. Moreover, as an 
intermediate step between the NJL and the a models 
a quark-meson model is obtained. Similarly, in the ap- 
proach of Cahill and Roberts [l3 . [Tsl it was shown by 
using bilocal auxiliary fields (along the same line of the 
Hubbard-Stratonovich transformation [l6j ) how to inte- 
grate out quark and gluon degrees of freedom to obtain 
a purely mesonic Lagrangian. While these calculations 
were only performed at one loop order, it was an inter- 
esting approach to connect effective models, say the a 
model, directly to QCD. In addition, there are also lat- 
tice QCD approaches for effective meson theories [iTj . 
which deliver qualitatively similar results. 

For the purpose of this paper it is necessary to chose 
a slightly different way, which consists of integrating out 
from the QCD Lagrangian the quarks only and then ob- 
tain a gluon-meson (fully bosonic) theory. From symme- 
try principles (colour gauge invariance and chiral sym- 
metry) we expect at leading order the following tree-level 
coupling between the gluons, the pions and their chiral 
partner, the scalar a meson: 

*-'gluon-meson ^ 

(G"'^''G;:j(a2+7f2 + ...) , (1) 

where C^^ is the gluonic field tensor and dots refer to 
other mesonic degrees, such as resonances with open and 
hidden strange-quark content (such as the kaons), vector 
resonances and, eventually, non-quarkonium resonances. 

The interaction Lagrangian ([T|) allows to study the 
transition of two gluons into a couple of mesons. If the 
latter are not stable, they further decay into pions or 
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kaons. For the mesonic sector we apply the most com- 
prehensive cr model, including not only scalar and pseu- 
doscalar mesons, but also vector and axial- vector mesons 
[igj . By assuming a thermal bath for gluons at the 
hadronization point, we can then evaluate how many pi- 
ous (and kaons) are obtained per gluon as function of the 
temperature T of the bath. To account for a possible fi- 
nite scale with an energy dimension, in the gluon sector 
at finite T, we allow for a finite effective gluon mass in 
the Boltzmann distribution. 

The paper is organized as follows: in Sec. II we derive 
our effective Lagrangian, in Sec. Ill we present the an- 
alytical and numerical results, and in Sec. IV we derive 
our conclusions. 



II. THE GLUON-MESON INTERACTION 

In this section we present the formal steps necessary 
to obtain a gluon-meson theory. We start from the eu- 
clidean QCD generating functional 



Zi'n,^]^ J VqVqVAVQVuj 



(2) 



X e 



- CQCD+CGF+CppG-riq-gV 



where 



Cqcd = q [l^D^ + m.) q + -G^G^^^ , (3) 

A" 

fj and 77 are the fcrmion sources, Cqf and Cfpg represent 
the gauge-fixing and the Faddcev-Popov ghost terms. In 
a suitable gauge [13, [13: and with no fermion sources, 
the generating functional of QCD can be written as 

Z = I VqVqVAe- /. A^^d-^CK , 

(4) 

D^iix -V)= j VCjVuVAAl{x)Al{y) (5) 



where 



is the exact gluon propagator which contains all gluon 
sclf-intcractions and gluon-ghost interactions but ex- 
cludes quark loops. 

To introduce the mesons we multiply Eq. (j4]) by 

1 = -^y" 2?$Pi?,,2?L^e-^-^— , (6) 



where 



: Tr [a^$ta^$] _ ^^Tr [$^$] (7) 
lTr[i?^, + L^,]+^Tr[i?^ + L^] 



is the quadratic part of the globally invariant linear a 
model Lagrangian with U{Nf)ii x U{Nf)L symmetry: 
and are Nf x Nf Hermitian matrices for the 
(pseudo) scalar, vectorial right-handed and left-handed 
mesonic degrees of freedom [l9|: 



(8) 



tj. 1 



Ru 



Vf,-A^, 
- d" R" , L""" 



In the previous expressions arc the N'j generators of 
U{Nf), S"' = %/2gt°g are the scalar, P'' = V^qi^/^fq are 
the pseudoscalar degrees of freedom, = y^qj^^f^q, 
A'^ — y/2q^^'^^t°'q are the vector and axial-vector micro- 
scopic quark currents. For instance, in the case Nf = 2 
the fields are given by 



$ = (cr + iriN) t° + (ao + iTr)t 



(9) 



where t° = 1/2, t' = Ti/2, and t, arc the Pauli matrices. 
Analogously, R^ = {uj^ — /C) + {p^ — a^)t represents 
the vector and L'" = (w'^ + /f) + (p^ + a^)f the axial 
vector degrees of freedom. The extension to the case 
A^y 3 is straightforward [20j . 

Clearly, upon Gaussian integration, Eq. © gives a 
constant factor and thus does not change the path inte- 
gral in the generating functional Eq. (j?]). 

To couple the mesons to fermions, the mesonic fields 
can be added as parallel transports in the fermion matrix 
similar to a mass term in the case of the scalar sigma 
and to a chiral transport in the case of the pion. The 
corresponding minimal coupling is given by 

j^D^ ^ D ^ j^D^-ci<P~C2{j^V'' +-f^j^A^) , (10) 

where ci and C2 are free parameters, and is the axial 
current. Finally, performing the Grassmann integration 
over the fermion fields, we obtain a purely bosonic theory 



Z = I V^VL^VRf^VA 



det [d] . 



The coupling of the mesonic and gluonic degrees of free- 
dom resides in the fermion determinant det [l)] , which 
cannot be computed analytically. However, by requir- 
ing local colour gauge invariance and chiral symmetry 
and restricting to lowest dimensionality of the interac- 
tion Lagrangian, we are left to the following Lagrangian, 



.'gluon-mcson 



(11) 

where a and b are couplings with dimension Energy ^ 
and describe the transition from two gluons to two 
mesons. To obtain the relation between the parameters 
a, b and the parameters ci, C2 introduced in Eq. (|10p 
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TABLE I. The parameters for each meson pair initially pro- 
duced and then decaying into pions and kaons are the mass, 
degeneracy, number of pions produced by the meson, number 
of kaons produced by the meson, the family factor. For me- 
son octets and their respective chiral partners we assume the 
same family factor. 



it would be necessary to evaluate the fermion determi- 
nant det [D] analytically. Although it is not possible to 
compute the fermion determinant exactly, it is natural 
to expect that a b. In the following we will work with 
the simplified assumption a = b. Note also that the La- 
grangian makes use of the so-called flavour blindness 
of the gluon fields. 

In accordance with Parganlija et al. [l^, the scalar 
partners of the pions are identified with the scalar res- 
onances above 1 GeV. However, for completeness, one 
should also include the scalars below 1 GeV, which ac- 
cording to many recent and less recent studies are non- 
quarkonium states. In fact, these states can be enhance- 
ments in the two-pseudoscalar channel or tetraquark 
states, see Refs. [2] and references therein for a more 
detailed discussion of this point. Here they are coupled 
to the gluons with an independent coupling c. We shall 
here test two choices: c = a, i.e. with equal strength as 
the other channels, and c = 0, where the light scalars are 
switched off. 

The explicit evaluation of the traces delivers 

•^gluon-mcson = aC^'^^ G°^j^{'k'^ + • • • + Oq + • • • ) 

+ bC^^^^^Gl^iPs -p' + ---+au-a^' + ■■■) 
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(12) 
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FIG. 1. Number of (top) pions Ilw{Epair) and of (bottom) 
kaons \lk{Epair) produced per gluon as a function of the gluon 
energy in the centre of mass of the gluon pair. Here a = b, 
c = are used. 



where dots refer to further quadratic mesonic interactions 
listed in Table H] 

Even if higher dimensionality terms in the gluon- meson 
Lagrangian exist, with more derivative or field terms, the 
reaction from two gluons to two mesons should domi- 
nate during the freeze-out at the boundary of the gluon 
plasma. Thus we employ this Lagrangian (fT2|) to evaluate 
the number of produced pions and kaons per gluon. 



III. ANALYTICAL AND NUMERICAL 
RESULTS 



In our framework an effective gluon mass Alg can be 
introduced. Notice the existence of a possible effective 
gluon mass 21|, or pole in a gluon propagator, already 
at T = 0, has been under debate in QCD. While gauge 
invariance in a perturbative approach rules out a gluon 
mass and the existing lattice QGD glueballs suggest that 
the gluon propagator is transverse, nevertheless a run- 
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ning gluon mass does not contradict gauge invariance 
[23 | , and there is also evidence both from Landau Gauge 
lattice QCD and from the glueball spectrum for a 
finite pole in the gluon propagator. Most results, sev- 
eral of them from lattice QCD calculations, point to a 
T = effective gluon mass or to other possible scales, 
say an effective dual gluon mass [l^l, in the range [0.5, 
1.0] GcV. A possible gluon mass at finite T is also start- 
ing to be investigated in lattice QCD pa . Moreover in 
the BAMPS set-up a finite Debye mass Q for the gluon 
is also considered, 



Mr, 



24 



(13) 



which, say at T = Tc = 0.158 GeV and as — 0.3 leads to 
a gluon mass of Mg ~ 0.239 GeV. 

For the following purposes we evaluate the Lorcntz- 
invariant Mandclstam variable s for a system of two glu- 
ons with four-momenta 



Pi 



A/2, 



1,2 



(14) 



The Mandelstam variable s leads to the centre of mass 
energy Epair of the gluon pair, 

s = (pi +P2)^ 
2 



(15) 
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Due to the Lagrangian in Eq. (fT2|) . the two gluons con- 
vert into a meson-pair. Considering that only pions and 
kaons are regarded as stable, we must also take into ac- 
count that the other resonances decay subsequently into 
kaons and pions. For instance, each a meson decays into 
a two-pion pair, therefore the aa channel results into a fi- 
nal 47r state. Similarly, a pair of p mesons decays also into 
four pions. On the contrary, the vector state (j) decays 
predominantly into kaons. In Table |I] these conversion 
factors, expressed as n^'^^ and •nS'^\ are listed for all the 
relevant mesons used for the calculations. Moreover, we 

also include the usual phase space factor v — Af|, 
where is the mass of the fc-th meson pair, and the 
corresponding degeneracy spin-isospin factor f/^ . Finally, 
we should also take into account the relative strength, 
which is equal to in the (pseudo)scalar channel, to 6^ 
in the (axial-) vector channel, and to in what concerns 
the scalar channel below 1 GeV. Putting all together, the 
number of pions and kaons per gluon is calculated as 



Mi0{Epa^r - ^Mk)gkfk^n\ 



(16) 

where i = tt, K. 

The functions TlT^{Epair) and Hk (Epair) as a function 
of the gluon energy and as a function of the Boltzmann 
temperature arc both depicted in Fig. [T] 



FIG. 2. Number of (top) pions and (bottom) kaons produced 
per gluon as function of the temperature T for different gluon 
masses. 



We are interested in comparing our results with the 
parametrization of the BAMPS set-up since, like our 
framework, BAMPS also consists of gluons, decaying into 
pions. In BAMPS, 1.5 to 2.0 pions are produced per 
gluon pair. This production takes place in different con- 
ditions than ours, at non-chemical equilibrium with local 
and dynamical many-gluon simulations. Nevertheless, if 
we compare with our approach, we notice we can eas- 
ily obtain Nj^ ^ 1, but the increase of Nt^ above unity 
can only be achieved at the price of including a sizeable 
effective gluon mass. This confirms the importance of 
including in our framework an effective gluon mass, sim- 
ulating a finite non-perturbative scale characteristic of 
the gluon plasma. 

We now evaluate the emitted number of pions N-^ (T) 
and kaons Nk{T) per gluon as a function of the tempera- 
ture T . We denote the multiplicities Ni{T) by employing 
a Boltzmann distribution of each gluon, thus leading to 



dpidp2fB{plT)fB{plT)n,{Epair) 



(17) 



5 




100 200 300 400 500 600 700 800 900 1000 
M„ [MeV] 



FIG. 3. We show (solid line) our result and (band) the BNL, 
CERN and FNAL data for the kaon to pion multiplicity ratio 
as a function of a possible gluon mass, at T = Tc ~ 158 
MeV. The measured ratio Nk/N^, £ [0.13, 0.22] is realized for 
a gluon mass Mg G [0.25, 0.35] GeV. 



where /s(p2,T) = 7Ve"VP''+"9'"°"/^ is the normahzed 
Boltzmann distribution. The integral in Eq. () can 
be simplified to a three-dimensional integration, and we 
compute it with a numerically accurate cH — h code. In 
Fig. [2] the functions N^{T) and Nk{T) are plotted for 
different values of the gluon mass Mg and for both choices 
a = 6, c = and a = b = c. For the intermediate value 
Mg = 400 MeV we have roughly one pion per each gluon 
for each T, while Nk{T) is a rapidly increasing function 
with T. 

In Fig. 13] we present the ratio Nk{T)/N^{T) for the 
temperature oiT ^ 0.158 GeV as a function of the 
gluon mass. The critical temperature Tc for the confine- 
ment and chiral crossover was measured in lattice QCD to 
be in the range € [0.145, 0.165] GeV [H-il. This was 
achieved in very precise full lattice QCD simulations with 
dynamical fcrmions. This temperature range is consis- 
tent with the freeze-out temperature of the quark-gluon 
plasma measured in Heavy ion collisions. The freeze-out 
temperature in heavy ion collisions can be determined 
from the inverse slope of the hadronic species multiplic- 
ity as a function of the transverse momentum. Recent 
analysis of heavy ion collisions indicate that the freeze- 
out temperature is in the range of [0.150, 0.170] GeV with 
results between 0.150 GeV and 0.160 GeV [3011 and results 
between 0.160 GeV and 0.170 GeV [MI,[32rBoth ranges 
are compatible, and we consider in our computations the 
mean value of T = 0.158 GcV. 

In Fig. [3] we compare the ratio with the experimen- 
tal data (A^A:/^ir)oxp € [1.13,1.22] measured by the 
PHENIX, STAR, BRAHMS, E866 and NA49 collabora- 
tions and extrapolated by the UrQMD 2.0, UrQMD 2.1 
and HSD transport approaches [H,!!^ for the ratio of the 



pion and kaon multiplicity in the most central collisions. 

Our results point to a solution corresponding to a pos- 
sible effective gluon mass, in a range of Mg S [0.28, 0.37] 
GeV. Remotely, a second less likely mass of circa 0.8 GeV 
may be possible. We notice that the solution points to 
a gluon mas at Tc of the order of circa 0.4 of the gluon 
effective mass of 0.5 to 1.0 GcV at T = resulting from 
different gluon calculations. The solution for the gluon 
mass is also consistent with the Debye mass Q of the 
gluon at finite T. 

Finally we test the robustness of our results checking 
the parameter dependence of the gluon-meson model. In 
Fig. U we compare the ratio Nk/Nt^ for different cou- 
plings (c = with c — 1), different temperatures (T ~ 
0.145 GeV with T = 0.170 GeV) and different statisti- 
cal distributions (Boltzmann and Bose-Einstein) . All the 
tests we performed with plausible changes of our param- 
eters suggest our results are robust. 

The small effect of changing the coupling of gluons to 
mesons is quite relevant for our results, since the meson 
properties are not yet established. Strong coupled chan- 
nel effects, tetraquarks, and giueballs have been shown 
to affect the meson spectrum and the mesonic couplings. 
Here we utilize the meson properties listed in the Parti- 
cle Data Group [131, but other meson data would yield 
similar results. We utilize the sigma model for the meson 
production, but other hadronic models would again yield 
similar results. 

The robustness of our results occurs because a gluon in 
the freeze-out of the plasma has a low energy (T ~ 160 
MeV and an effective mass of Mg ~ 320 MeV), much 
lower than the energy of a gluon in any glucball typical 
of lattice QCD simulations or of constituent gluon model 
estimations. Thus our results do not really depend on 
the details of the meson spectrum and of the meson cou- 
plings above that low energy, and escape the problem of 
understanding higher energy reactions such as the glue- 
ball decays. 



IV. CONCLUSIONS 

In this work we develop an effective Lagrangian which 
connects gluons to mesons. We then utilize this approach 
to calculate the emitted number of pions and kaons per 
gluon out of a gluon gas at temperature T. We assume 
the gluons are in a thermal bath at T — Tc whereas the 
mesons are produced at vanishing temperature. 

The fact that our effective Lagrangian consists of 
bosonic degrees of freedom only might be also interest- 
ing for lattice QCD applications. In fact, the Grassmann 
variable nature of the quarks makes them technical!; 
much more difficult to address (see for instance Ref. 
and refs. therein) than the bosonic gluons and mesons. 

Our approach represents an attempt to link directly 
and in an understandable way the gluon-dominated 
physics of the quark gluon plasma, as suggested by the 
colour glass condensate and BAMPS approaches, to the 
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FIG. 4. The parameter dependence in the gluon-meson model. We compare the ratio Nk/Nt^ for different coupUngs (left, c = 
with c = 1), different temperatures (centre, T = 0.145 GeV with T — 0.170 GeV) and different statistical distributions (right, 
Boltzmann and Bose-Einstein) . 



light hadrons in the final stage. 

Further developments of our approach could include 
meson mass modifications in the medium, effects of 
the freeze-out boundary, different meson-gluon couplings, 
missing resonances such as the tensor mesons, final state 
interactions among mesons, and, last but not least, glue- 
ball fields, which directly couple to gluonic fields. Never- 
theless, all the tests we performed with plausible changes 
of our parameters suggest that our results arc stable. 
Thus we regard our effective Lagrangian as a pilot study 
toward a better understanding of the rich physics of QCD 
and of Heavy Ion Collisions. 

Interestingly, to reproduce the experimental results of 
the PHENIX, STAR, BRAHMS, E866 and NA49 coUab- 
orations for the pion and kaon multiplicities, we have to 
include a finite scale for the gluon energy at T = T^. Our 
results point to a possible effective gluon mass, in a range 
of Mg between 0.28 and 0.37 GeV if we consider the un- 
certainty on the rapidity ratio, or of Mg G [0.25, 0.39] if 
we also consider the uncertainty on the freeze-out tem- 
perature. We notice that this solution points to a gluon 
mass at Tc of the order of circa 0.4 and 1.0 of the gluon 
effective mass of 0.5 to 1.0 GeV at T = resulting from 
different gluon calculations. An Mg € [0.25,0.39] is also 



close to the Debye gluon mass at T = Tc. Notice this 
solution corresponds to an absolute pion multiplicity of 
circa 0.9 pious per gluon. 

If the gluon mass is related to confinement, say as in 
superconductors, our result is consistent with the QCD 
order parameters at Tc. Both for a first order phase tran- 
sition where the order parameter is discontinuous (found 
in quenched lattice QCD [s^) and in a crossover where 
the order parameter remains finite (found in lattice QCD 
with dynamical fermions [26j). the scale of confinement 
should not simply vanish at T = Tc. Here we present 
an evidence, based in the experimental Heavy Ion data 
and in our simple and robust gluon-meson model that a 
finite scale of 0.25 to 0.39 GeV exists in the gluon sector 
of QCD at T = Tc. 
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